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Abstract. This paper is devoted to the study of a nonhnear coupled system 
of two Korteweg-de Vries equations in a periodic domain under the effect of 
an internal damping term. The system was introduced Gear and Grimshaw to 
model the interactions of two-dimensional, long, internal gravity waves prop- 
agation in a stratified fluid. Designing a time-varying feedback law and using 
a Lyapunov approach we establish the exponential stability of the solutions in 
Sobolev spaces of any positive integral order. 



1. Introduction 

The goal of this paper is to investigate the decay properties of the initial- value 
problem 

u' + uux + Uxxx + as,Vxxx + aivvx + a2{uv)x + k{u - [u]) = 0, 
hiv' + rvx + vvx + Vxxx + haaUxxx + b2a2UUx 

+b2ai{uv)x + k{v - [v]) = 0, 

u(0, x) = 0(a;), 
^(O, x) = ip{x) 

with periodic boundary conditions. In (11.11) . r, ai, 02, 03, 6i, &2, are given real 
constants with 61, 62, fc > 0, u{t, x), v(t, x) are real- valued functions of the time and 
space variables t > and < x < 1, the subscript x and the prime indicate the 
partial differentiation with respect to x and t, respectively, and [/] denotes the 
mean value of / defined by 

[/] := /' fix) dx. 



(1.1) 



^0 

When k — 0, system was proposed by Gear and Grimshaw [5] as a model to describe 
strong interactions of two long internal gravity waves in a stratified fluid, where the 
two waves are assumed to correspond to different modes of the linearized equations 
of motion. It has the structure of a pair of KdV equations with both linear and 
nonlinear coupling terms and has been object of intensive research in recent years. 
In what concerns the stabilization problems, most of the works have been focused 
on a bounded interval with a localized internal damping (see, for instance, [14^ and 
the references therein). In particular, we also refer to [1] for an extensive discussion 
on the physical relevance of the system and to [3l EJ [3 [6l [7] for the results used in 
this paper. 

We can (formally) check that the total energy 

1 

E = 7: hu'^ + biv'^ dx 
2 Jo 
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associated with the model satisfies the inequality 

E'^~k[ b2{u - [u]f + {v ~ [v])^ dx < 
Jo 

in (0, oo), so that the energy in nonincreasing. Therefore, the following basic ques- 
tions arise: are the solutions asymptotically stable for t sufficiently large? And if 
yes, is it possible to find a rate of decay? The aim of this paper is to answer these 
questions. 

More precisely, we prove that for any fixed integer s > 3, the solutions are 
exponentially stable in the Sobolev spaces 

H;{0, 1) {u e H%0, 1) : d^u{0) = d'^u{l), n = 0,...,s} 

with periodic boundary conditions. This extends an earlier theorem of Davila in 
[6] for s < 2. 

Before stating the stabilization result mentioned above, we first need to ensure 
the well posedness of the system. This was addressed by Davila in [3] (see also [1]) 
under the following conditions on the coefficients: 

a§62 < 1 and r = 
Q 2) ^20103 - bias + bia2 - 02 = 

fcifli — fli — 610203 + 03 = 
6ia| + 6201 — biai — a2 — 0. 

Indeed, under conditions (|1.2p . Davila and Chaves [7] derived some conserva- 
tion laws for the solutions of (ll.ip . Combined with an approach introduced in 
[H [T7] , these conservation laws allow them to establish the global well-posedness in 
Hp{0, 1), for any s > 0. Moreover, the authors also give a simpler derivation of the 
conservation laws discovered by Gear and Grimshaw, and Bona et al 1 . We also 
observe that these conservation properties were obtained employing the techniques 
developed in [T3] for the single KdV equation; see also [T^ . 
The well-posedness result reads as follows: 

Theorem 1.1. Assume that condition p.2p holds. If (f),4' G ^^{0,1) for some 
integer s > 3, then the system (|l.ip has a unique solution satisfying 

u, V e c([o, 00); h;{o, 1)) n ^^([o, ^)-h;-\o, i)). 

Moreover, the map ((/i, "0) ^ {'^1'^) continuous from (i7p(0, 1))^ into 

(c([o, 00); h;{o, 1)) n ^^([o, ^);h;-^o, . 

For fc = 0, the analogous theorem on the whole real line — 00 < a; < c» was 
proved Bona et al. [T], for all s > 1. 

With the global well-posedness result in hand, we can focus on the stabilization 
problem. For simplicity of notation we consider only the case 

(1.3) 6i = 62==l. 
Then the conditions (|1.2p take the simplified form 

(1.4) r = 0, a\ + a\ = ai ~\- 02, |a3| < 1, and (ai — 1)03 = [02 — 1)03 — 0. 

Hence either a^ = Q and a\ + a\ = ai + 02, ot: < \as\ < 1 and ai = 02 = 1. 
We prove the following theorem: 

Theorem 1.2. Assume (|1.3p and (|1.4p . If <f),^ G Hp{0, 1) for some integer s > 3, 
then the solution of (|l.ip satisfies the estimate 

Mt) - ["(<)] IIh=(o,i) + Mt) ~ bW]llH=(o,i) = o ' * ^ °o 

for each k' < k. 
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An analogous theorem was proved in \10\ for the usual KdV equation by using the 
infinite family of conservation laws for this equation. Such conservations lead to the 
construction of a suitable Lyapunov function that gives the exponential decay of the 
solutions. Here, we follow the same approach making use of the results established 
by Davila and Chavez [7]. They proved that under the assumptions (|1.2p system 
(jl.ip also has an infinite family of conservation laws, and they conjectured the 
above theorem for this case. At this point we observe that some computations are 
simplified if we change u, v, cf) and tp to u — [u], v — [v], cf) — [0] and tp — Then, 
the new unknown functions u and v satisfy the same system with ku and kv 
instead oik{u~ [u] ) and k{v—[v\). Hence we consider the solutions of the simplified 
system 




u' + uux + Uxxx + a'iVxxx + aivvx + a2{uv)x + fcu = 0, 
v' + vvx + Vxxx + asUooxx + a2UUx + ai{uv)x + few = 0, 
u(0, x) — (f>{x), 
u(0, x) = ip{x) 



with periodic boundary conditions, corresponding to initial data (j), ip with zero 
mean values. 

In order to obtain the result, we prove a number of identities and estimates for 
the solutions of (II. ip . In view of Theorem 1 1.11 it suffices to establish these estimates 
for smooth solutions, i.e., to solutions corresponding to C°° initial data 0,?/' with 
periodic boundary conditions. For such solutions all formal manipulations in the 
sequel will be justified. 

Finally, we also observe that a similar result was obtained in [\\^ for the scalar 
KdV equation on a periodic domain. The authors study the model from a control 
point of view with a forcing term / supported in a given open set of the domain. 
It is shown that the system is globally exactly controllable and globally exponen- 
tially stable. The stabilization is established with the aid of certain properties of 
propagation of compactness and regularity in Bourgain spaces for the solutions of 
the corresponding linear system. We also refer to |llj for a quite complete review 
on the subject. 

The paper is organized as follows. In Section 2 introduce the basic notations 
and we prove some technical lemmas. Sections 3 to 6 are devoted to the proof of 
the exponential decay in H!l, for s = 0, 1, 2 and s > 3, respectively. 



2. Some technical lemmas 

In the sequel all integrals are taken over the interval (0, 1) so we omit the inte- 
gration limits. 

As explained in the introduction, all integrations by parts will be done for smooth 
periodic functions. Therefore, we will regularly use the simplified formulas 

j 1x9 dx ^ ~ j fgx dx and J f" fx dx = (n==0, 1,...) 

without further explanation, and we will also use the simplified notation 

d" f 

As an example of the application of these rules we show that the mean- values of of 
the solutions are conserved: 



Lemma 2.1. The mean-values [u] and [v] of the solutions of (jl.ip do not depend 
on t. 
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Proof. We have 

= - f uux + Uxxx + a^iVxxx + aivvx + a2{uv)x + A;(m - [u]) dx 
+ Uxx + azVxx + "1 Y + a2'«wj + - M) dx 

= —k J {u — [u]) dx 
= 

and 

H' = vxxx + a.uxxx + a^uux + a, + " H) dx 

= —A; J ^ b]) 
= 

by a straightforward computation. □ 
Theorem 11.21 will thus follow from the following proposition: 



Proposition 2.2. Under the assumptions of Theorem ] l.SA the smooth solutions of 
3)) satisfy the identity 



(2.1) J u{t)^ +v{tf dx = e-^''^ J <j)^ dx, t>0, 
and the estimates 

(2.2) e^*^'* J {d^u{t)f + {dy{t)f dx ^ as t ^ oo 
for all positive integers n and for all k' < k. 

Remark. For n = I the proposition and its proof remain valid under the weaker 
assumption that jasj < 1. We can also add the term rvx to the equation by changing 
g to g — rv^ in Lemma |4. II 

Proposition 12.21 is proved by using the Lyapunov method. More precisely, we 
shall use the following lemma: 

Lemma 2.3. Let f : (0, cxd) ^ R &e a nonnegative function, and write hi « h2 if 
hi — h2 = o(f) as t oo. 

// there exists a function g : (0, oo) —5- M such that g « 0, f + g is continuously 
differentiable, and (/ + g)' ~ —"ikf for some positive number k, then 

e^^'^f{t) ^0 as t^oo 

for each k' < k. 

Proof. Fix k" > such that k' < k" < k, and then fix e > such that 

1 -£ _ k^ 

1 + £ ^ k ' 

Finally, choose a suflaciently large <' > such that 

(l-e)/(i)<(/ + 5)(t)<(l + £)/W 

and 

2k{l - e)f{t) < -if + gYit) < 2fc(l + e)f{t) 
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for all t > t' . Then for t > t' we have 



-(/ + gYit) > 2fc(l - e)f{t) > 2k^{f + g){t) = 2k" {f + g){t), 



whence 



It follows that 



for all t > t' , and hence 



|(e--(/ + <0. 
5''"*(/ + 5)(t)<e^^"*'(/ + 5)(t') 



< e^^'^f{t) < 



„2k"t' 



"''{f + 9W) -m"-k')t 



1 -e 



for all t > t'. We conclude by observing that e ■^'^'^ — as i — cx). 



□ 



For the proof of the next result, we shall use the Holder and Poincare-Wirtinger 
inequalities in the following form. The second estimate will be used only for func- 
tions with mean value zero: [u] =0. 

Lemma 2.4. Ifp,q e [0,oo), then 

(2.3) < for all m6L'?(0,1) and 1 < p < q < oo; 

(2.4) ||w - Mllp < llwxilg for all ueH^{0,l) and l<p,q<oo. 

We shall frequently use Lemma [2.31 together with the following result: 

Lemma 2.5. Let n > I and let armPm, fn — 0, ...,n be nonnegative integers 
satisfying the two conditions 

2(a„ + /3„) + a„_i +/?„_! <4 



and 



Then 



m=0 



m=0 

If, moreover, d> 3 and 



< 



+ dx 



+^'n-l dx 0, 



then it follows that 

as t ^ oo. 
Proof. Setting 

we have 



J n dx = o(^j ul + vl dx^ 



a/u^ + and ^m-=am + l3m, m = 0,...,n 



/n 
™ — n 



< 



n 



dx. 



We are going to majorize the right side by using the Holder and Poincare-Wirtinger 
inequalities (|2.3p - (|2.4p . We distinguish five cases according to the value of 7„ +7,1-1: 
since 27„ + 7,1-1 < 4 by our assumption, 7„ + 7„_i < 4. 
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If 7n + 7ri-i = 0, then we have 



/n 
n 



< 



ni 

m=0 



\Zm\\l^ < \\Zn\\2 ll^n-lUa' 



If 7n + 7n-i = li then 



/n 
n 



m— 

If In + 7"-i = 2, then 



n-2 



< ll^"|li n ^ Ik"|l2 I|2«-1||2 

m— 



d-2 



/n 
n 



n|l2 



ji-2 

ni 

m=0 



< \\Zn\\l \\Zn-l\\i ^ 



If 7n + 7ri-i = 3, then we have necessarily 7n = 1 and 7„_i ~ 2, so that 



/n 
n 

m=0 



kmlli^ < ll^nlls I|Z„-1||2 



d-2 



< ii^«ii2pn-iiLp«-iii2 n 

m=0 

Finally, if 7„ +7,1-1 = 4, then we have necessarily 7^ = and 7„_i = 4, so that 

ri-2 



n '^^ 

m— 



< Ikn-lIlL Il^»-l|l2 n - Il^"ll2 Il^"-l|l2 ^- O 



m— 



3. Proof of Proposition 12.21 for n = 
Our proof is based on the following identity: 
Lemma 3.1. The following identity holds for all n — 0,1, ... : 



(3.1) 



ul + vl dx 



-2k I + dx 



2 j Un{uiu)n+Vn{viv)n dx 
2ai J Un{vVi)n + Vn{uv)n+1 dx 
202 / Vn(uUi)n + Un(uv)n+1 dx. 



Proof. We have 



ul + vl dx 



2UnU'n + 2Vnv'j^ dx 

-~2un{uui + 1*3 + a^vs + aivvi + a2{uv)i + ku)n dx 
+ / -'2i;„(oti +V3 + asus + a2UUi + ai{uv)i + kv)n dx. 



This yields the stated identity because 

-2u„u„+3 - 2w„u„+3 dx ^ I 2u + 2Vn+lVn+2 dx 



and 



as 



+ dx^O 

j -2unVn+3 " 2w„u„+3 dx ^ j -2u„w„+3 + 2vn+zUn dx = 0. □ 
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Proof of the proposition for n = 0. In this case the last three integrals of the iden- 
tity p.ip vanish because 

uuiu + vviv dx ~ — 1 {u^ + v^)i dx = 0, 



uvvi + v{uv)i dx = j {uvv)i dx = 
and 

J vuui + u{uv)i dx = j {vuu)i dx = 0. □ 

Proceeding by induction on n, let n > 1 and assume that the estimates 
(3.2) ul^+vl^ dx = o (e"^'''*) as t ^ oo 



hold for all integers m — 0, . . . ,n — 1 and for all k' < k. For n = 1 this follows from 
the stronger identity (|2.ip . 

4. Proof of Proposition 12.21 for n = 1 

For the proof of the case n = 1 we shall use an identity suggested by a conser- 
vation law discovered by Bona et al. [I]. 

Lemma 4.1. Setting 

r 

f := ul + vf + 2a3UiVi dx 

and 

5 — - / {u'^ + v^) + 3{aiuv^ + a2y?v) dx, 



3 „ 

we have the following identity: 

(4.1) (f + g)' = -2kf ~ikg. 

Proof. The equality (|4.ip will follow by combining the following four identities: 



(4.2) 



j Ui+vf dx^ = —2k J Ui+Vi dx 

— y ul + vf dx 

— 3ai J uiv'l dx 

— 3a2 J ulvi dx; 



(4.3) ( / uiVi dx ) = —2k I uivi dx + j uuiV2 + vviU2 dx 



Ol 

2 



J 2v2Uiu -t- 3wiu^ -|- v'l dx 

0,2 j 2^7 

— / 2u2Viv + 6UiVi + Ui dx; 
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(4.4) [ I + dx] = -3k j + dx ~ i j u\ + vl dx 

— ai j iu^vvi + 2v^ui dx 
r 

— 02 / Sv^uui + 2u'^vi dx. 
+ 6a3 / UU1V2 + VV1U2 dx; 



(4.5) y^J aiuv + v dx j — — 3fc J aiuv + v dx 

+ ai J g'^^'^'^'i + u^vvi — Sv^ui dx 

/2 
-M'^wi + v^uui — 3uiVi dx 

— aifls J 2v2Uiu + 3viUi + vf dx 

— 0203 y 2u2V\v + 3u]v\ + dx. 

Proof of (|4.2p . We transform the identity p.ip for n = 1 as follows. We have 

ui{uiu)i + vi{viv)i dx — J U2U1U + u\ + V2V1V + vf dx 

f 1 1 

= ul + vf + -{uDiu + -{vf)iv dx 

= \ ul + vl dx, 



ui(vvi)i + vi{uv)2 dx = J uivf + U1VV2 — V2{uv)i dx 

= J uiv\ — V2UV1 dx 

uivl - ]^u{vl)i dx 
2 



uiv-y dx 



and by symmetry 



vi[uui)i + ui{uv)2 dx — — I u\vi dx. 



Using them (P?T|) implies 
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Proof of (1131). We have 

J uivi dx^ = J u[vi + uiv'i dx 

= J —[uui + Uz + a^vj, + aivvi + a2{uv)i + ku)iVi dx 

+ j —ui{vvi + -^3 + 031*3 + a2UUi + ai(u-i;)i + kv)i dx 
= —2k J uiVi dx + J {uui + U3)v2 + {wi + V3)u2 dx 

— ai J (wi)ifi + ui{uv)2 dx 

— 02 y {uv)2Vi + ui{uui)i dx 

— aa y + U4U1 dx 

= —2k J uivi dx + J UU1V2 + VV1U2 dx 
+ ai y VV1V2 + U2{uv)i dx 
+ a2 {uv)iV2 + U2UU1 dx 



because 



U3V2 + V3U2 dx = / U3V2 - V2U3 dx ^0 



and 



V4V1 + U4U1 dx = ~ I V3V2 + U3U2 dx — —- I {v2 + ^2)1 dx — 0. 



Since 



/ 



VV1V2 + U2{uv)i dx — j -v{vl)i + -{ul)iv + U2UV1 dx 



f 1 1 

/ ^2^1 ^ 2^^^^ ~ ^^^^ ~ U1UV2 dx 



2 

and by symmetry 



1 I 2 3 7 

- / 2v2Uiu + SviUi + Vi dx 



UU1U2 + V2{uv)i dx = / 2u2Viv + iuivl + u\ dx. 



2 

(|4.3p follows from the previous identity. 
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Proof of (|4.4p . We have 

dx ] — / 3u^u' dx 



— 3u^(uui + U3 + 03 + aivvi + a2{uv)i + ku) dx 
y — — (it^)-^ + 3u {u'l) ^ — iku^ dx — 3a3 j u^v^ dx 



— 3ai J u^wi dx — 3a2 j u^vi + -^{u^)iv dx 
= —3 j u\ + ku^ dx — 3ai j u^vvi dx — 2a2 J u^vi dx 
+ 603 / UU1V2 dx. 



We have an analogous identity for J dx by symmetry; adding the we get (j4.4p . 
Proof of (|T5)) . We have 



-2uv{uui + U3 + 03(73 + aivvi + a2(Mw)i + few) dx 

+ J —u^{vvi + W3 + 037/3 + a2UUi + ai{uv)i + kv) dx 
2u^uiv + 2u2{uv)i — u^vvi + 2v2uui dx — 3k J v}v dx 

— ai J 2uvvvi + u'^{uv)i dx 

— a2 J 2uv{uv)i + u^ui dx 

— 03 / 2uvv^ + M^U3 dx. 



Here 



J —2v}u\v dx = — - J {u'^)iv dx ~ — J u^vi, 

^ ' .2l„.2\ , 1 / /„,2^ „2 , /„,2„ 



-u vvi dx — —— u {v )i dx — — I {u )iv dx = j v uui dx, 



2u2{uv)i + 2v2uui dx ~ J {2u2Uiv + 2u2uvi) — {2viu1 + 2V1UU2) dx 
= J {ui)iv — 2viul dx 
= —3 J ulvi dx, 

J 2uvvvi + u^{uv)i dx — J —u{v^)i + u^vi + —{u'^)iv dx — — J u^vi — v^ui dx, 
2uv{uv)i + u'^ui dx — j ({uv)^ + '^^j '^•^ ~ 
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and 

2uvv3 + u^u^ dx ^ I ~'2{uiv + uvi)v2 ~ 2uuiU2 dx 

2{u2V + uiVi)vi — u{v\)\ — dx 
2{u2V + uivi)vi + uivf + ul dx 
2u2Viv + 3uivf + ul dx, 
so that 

u^v dx^ — J g'^^'^^i + v^uui — iu^vi dx — 3k J v?v dx 

ai j u'^vi — v^ui dx — j 2u2Viv + 3uivJ + dx. 



3 

By symmetry, we also have 



v^u dx] — I —v^ui + u^vvi — 3v\ui dx — 3k j v'^u dx 



~ -^0,2 j V ui — u vi dx — J 2v2Uiu + 3viu^ + Vi dx. 

Combining the last two identities (|4.5p follows (some terms annihilate each 
other). □ 

Proof of the proposition for n — 1. It suffices to show that the functions / and g of 
Lemma satisfy the conditions of Lemma l2.3l Since [03 1 < 1, we have / > 0. The 
other conditions follow from the already proven case n = and from the second 
part of Lemma 12.51 We conclude by applying the lemma and then by observing 
that 

u\ + vf dx < ; — - / u? + + 2a3UiVi dx. □ 

1 - Iflal J 



5. Proof of Proposition 12.21 for n = 2 
Lemma 5.1. Setting 
r 

/ := / 1*2 + w| + 2azU2V2 dx, 

g :— —— J {uiU + vfv) + ai{2uiviv + vfu) + a2{2uiViu + ufv) dx 

and 

h := -as y (1 ~ ai){2u3V2U + U2V2U1) + (1 - a2){2vzU2V + U2V2V1) dx, 
we have 

(5.1) (^f + gY ^ -2kf + h. 

Proof. The relationship (|5.ip will follow by combining the following relations: 
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(5.2) ( I ul+V2 dx] = -2k I u^ + v^ dx 



5 / U2U1 + V2V1 dx 



5ai J 2u2U2'i'i + V2U1 dx 
5a2 / 2U2V2U1 + u\vi dx; 



(5.3) U2V2 dx^ — —2k J U2V2 dx 



(5.4) / ulu + vlv dx 



U3V2U + V3U2V + 3U2V2{UI + vi) dx 

ai J —{ul + V2)vi + 2U2V2U1 — U3V2U dx 
/ ^{"^2 + ^"2)^1 + 27/2^2^1 - V3U2V dx; 



« —3 J u\ui + v^vi dx 

— 2a3 / u^V2U + W3U2W + 2u2V2{ui + vi) dx; 



(5.5) 2uiViv + v\u dx^ Ki —2> j 2u2U2Wi + W2W1 dx 

+ —3(^2 + V2)vi + 2U3V2U — 2U2V2U1 dx; 



(5.6) [ / 2uiViu + u\v dx ) « —3 j 2u2V2Ui + u\vi dx 

+ az i —3(^2 + V2)ui + 2V3U2V — 2U2W2W1 dx. 



Proof of (|5.2p . Wc transform the last three integrals of the identity p.ip in the 
following way: 

— 2 J U2{uiu)2 + V2{viv)2 dx — —2 J 3M2M1 + U2U3U + 3V2V1 + V2V3V dx 

= —2 J 3M2W1 + -(^2)1" + 3V2V1 + —{v2)iv dx 
r 

= — 5 / U2U1 + V2V1 dx, 
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2ai / U2(oti)2 + V2iuv)3 dx = -2ai / 3u2ViV2 + U2VV3 - V3{uv)2 dx 



= — 2ai j ?>U2ViV2 ~ 2V3U1V1 — V3UV2 dx 
= -2ai / 3w2WiW2 + 2v2{uivi)i - \u{v2) 



f 5 2 

= — 2ai / 5M2W1W2 + 2'"i^2 cf^; 



-5ai / 2W2W2W1 + v^ui dx, 



and by symmetry 

— 2a2 J V2{uui)2 + U2{uv)3 dx = — 5a2 y 2u2W2Ui + u^vi dx. 

Combining these identities with p.ip we obtain (|5.2p . 
Proof of lEM- We have 

y U2V2 dx^ — J U2V2 + U2v'2 dx 

{ulU + U3 + ku + aswa + aiViv + a2{uv)i)2V2 dx 

U2{viv + V3 + kv + 03^3 + a2Wiu + ai(uu)i) 
2k I U2V2 dx 

«3 j V5V2 + U2U5 dx - J U5V2 + U2V5 dx 



{uUi)2V2 +W2(wi)2 dx 
fli J {vVi)2V2 +U2{uv)3 dx 
a2 / {uv)3V2 +U2{uUi)2 dx. 



Here 



J V5V2 + U2U5 dx = - J V4V3 + M3U4 dx = - - y (^3 + ^3)1 = 0, 



U5V2 + U2V5 dx — U5V2 — U5V2 dx — 0, 



[uUi)2V2 + U2{vVi)2 dx 

= / SU1U2V2 + UV2U3 + VU2V3 + iViV2U2 dx, 
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ivVi)2V2 + U2{uv)^ dx 



= I 3^2 wi + v:}V2V + U3U2V + 3U2V1 + 3U2V2U1 + V3U2U dx 

3u|ui + -{v2)iv + -^{uDi^ + SujWi + 3U2V2U1 + V3U2U dx 
5 

-(m2 + 1^2)^1 + 3U2V2U1 + V3U2U dx 
5 

-(W2 + V2)vi + 3U2V2U1 — V2U3U ~ V2U2U1 dx 

5 

— (^2 + V2)vi + 2U2V2U1 — U3V2U dx. 



By symmetry, we also have 

{uUi)2U2 + V2{uv)3 dx = f ^(^2 + 1^2)"! + 2u2V2Vi - V3U2V dx. 



This proves 



Henceforth in aU computations we integrate by parts and we apply Lemma 
several times. 

Proof of (I5.4P . We have 

,2„, \ _ / o ' „, 1 „,2„,/ 



dx j — I 2uiu'iU + u^u' dx 
~u'{2u2U + uf) dx 

{2u2U + u\){uiu + U3 + fcu + aiviv + a2{uv)i + a^vj^) dx 
k I 2u2'u^ + dx 

uiu{2u2U + uf ) dx 
U3{2u2U + uf) dx 
+ ai J viv{2u2U + Ui) dx 
+ 0-2 J {uv)i{2u2U + u1) dx 
+ 0-3 V3{2u2U + uf) dx. 



Here all integrals are equivalent to zero by Lemma |2.5[ except those containing U3 
or V3. Since 



U3{2u2U + uf) dx = {U2) lU + U3U1 dx = — / u\ui 2u\ui dx = —i I u\ui dx 
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and 

J vz{'2.U2U + u\) dx = 2 j vj,U2U — V2U2U1 dx 



= 2 y —V2U3U — V2U2U1 — V2U2U1 dx 
= —2 / U3V2U + 2U2V2U1 dx, 



we conclude that 

dx^ ~ —3 J U2U1 dx — 2a3 J U3V2U + 2u2V2Ui dx. 

Adding this to the analogous relationship for / vfv dx we get (|5.4p . 
Proof of ([531) and ^^6^. We have 

uiviv dx^ — J u'lViv + uiv[v + uiviv' dx 
= j —u'{v2V + v'l) — v'u2V dx 

{v2V + vf){uiu + U3 + ku + aiviv + a2{uv)i + a^v^) dx 

+ J U2v{viv + v:i + kv + a2Uiu + ai{uv)i + a3U:i) dx 
V2VU3 + vfus + U2VV3 dx + 03 J {v2V + vl)v-:i + U2VU3 dx 

{u2V2)lV — U2{Vi)i dx + J {V2V + vf)v3 + U2VU3 dx 

= —3 J U2V2V1 dx + as J {v2V + vf)v3 + U2VU3 dx. 

Since 

-{vDiv - 2vlvi + -viv^)! dx 
— — 2^2 wi — —U2V1 dx 



= / ^2^2^! - 2^2^! dx, 



it follows that 



J 2uiViv dx^ ~ —6 J U2V2V1 dx — J (5t;| + ^2)^1 dx, 
and then by symmetry 

2uiViu dx j —6 / U2V2U1 dx — I {^u\ + V2)ui dx. 
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Next we have 

J ufv dx^ = J 2uiu'iV + Uiv' dx 

— (2u2V + 2uiVi)u' + u\v' dx 

{2u2V + 2uiVi){uiu + W3 + /cm + aiViv + a2{uv)i + a^v^) dx 

+ J -~u1{viv + W3 + kv + a2Uiu + ai{uv)i + 03^3) dx 
2U3U2V + 2uiViU3 — ulv3 dx 

+ 03 y {2u2V + 2miWi)u3 — u\u^ dx 
—u\vi — 2u2(uiWi)i + 2U1U2V2 dx 

+ 03 y {2u2V + 2uiVi)v3 — ufu3 dx 
= —3 / u^vi dx + / {2u2V + 2uiVi)v3 — ulu^ dx. 



Since 

{2u2V + 2uiVi)v3 — dx = J —2v2{u^v + 2u2Vi + U1V2) + 2u'^ui dx 

r 

= / — 2U.3U2W — — 2V2U1 + 2u\ui dx 

= 2 y U3U2W - U2V2V1 + {u\ - vl)ui dx, 

it follows that 

dx ) = —3 / U2V1 dx + 2a3 / V3U2V — U2V2V1 + (uj — V2)ui dx, 



and then by symmetry 

J v\u dx^ = —3 j V2U1 dx + 203 j U3V2U — U2V2U1 + — ^2)^1 
Combining the four relations we get (|5.5p and (|5.6p . □ 

Proof of the proposition for n = 2. We consider the functions f,g,hoi Lemma l5.ll 
If 03 = or if fli = 02 = 1, then h — Q. If [03 1 < 1, then 

+V^dx< j r / u\-\-V^ + 2a3UnVn dx . 



Since by Lemma 12.51 and the induction hypothesis / and g satisfy the assumptions 
of Lemma 12.31 we may conclude as in case n — 1 above. □ 

6. Proof of the proposition for n > 3 

We proceed by induction on n, so we assume that the proposition holds for 
smaller values of n. 
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By Lemma |3. II we have 
(6.1) (^J ul + vl dx^ = -2fc J ul + vl dx 



2 y Un{uxu)n+ Vn{vxv)n dx 
2ai J Un{vVi)n + Vn{uv)n+1 dx 
2a2 / Vn{uUi)n + Un{uv)n+1 dx. 



If we differentiate the products in the last three integrals by using Leibniz's rule and 
the binomial formula, we obtain a sum of three-term products. Using the inequality 
n > 3, it follows from Lemma |2 . 5 1 that all terms are equivalent to zero, except those 
containing the factor Un+i or w„+i. 

Indeed, the orders of differentiation of the three factors are n, j and n + 1 — j 
with 1 < j < n. Since the sum 2n + 1 of the differentiations satisfies the inequality 
2n + 1 < 2n + (n — 1), we have 

2 (a„ + /3„ ) + +Pn-i) <4, 

and Lemma 12.51 applies. 

Using again that 1 < n — 2, it follows that 

Un{uiu)n + Vn{viv)n dx M„U„+iU + VnVn+lV dx 



uf-^ui + v'^vi dx 



\ j {Un)lU+{vl)iV dx 
1 

^2 
0, 



Un{vVi)n + Vn{uv)n+1 dx ^ j M„W„+i + W„M„+iU + W„UW„+i dx 

UnVVn+1 - Un{v,iV)i + dx 
= J -UnVnVi - ^Uivl dx 

and by symmetry 

J U„(uUi)„ + Un{uv)n+1 dx W 0. 

Using these relations we infer from (|6.1|) that 

u^, + vl dx] « —2k I ut +v'i dx 



and we conclude as usual. 
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